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In this work, we present numerical resuhs which support 5*0(5) symmetry as a concept unifying 
superconductivity and antiferromagnetism in the high-temperature superconductors. Using exact 
cluster diagonalization, we verify that the low-energy states of the t — J model, a widely used 
microscopic model for the high Tc cuperates, form SO{5) symmetry multiplets. Our results show 
that the d— wave superconducting ground states away from half-filling are obtained from the higher 
spin states at half-filling through 5*0(5) rotations. 
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Introduction: Ever since the early days of quantum mechanics, group theoretical interpretation of spectroscopy 
revealed deep symmetry and profound unity of Nature. Atomic spectra can be fitted into irreducible representations 
(irreps) of SO {3), and the regular patterns which emerged from this classification offered fundamental understanding 
of the periodic table. After the discovery of a large number of hadrons, the "embarrassment of riches" is removed 
by the classification of hadronic spectra into irreps of SU{5) and this hidden regularity inspired the predictions of 
quarks, the fundamental building block of the universe. In our quest for understanding the fundamental design of 
Nature, the importance of symmetry can never be over-emphasized. 

In this article we report a different kind of spectroscopy, and its classification into a different kind of symmetry. The 
spectroscopy is performed on a computer, which numerically diagonalizes microscopic Hamiltonians, i. e. Hubbard 
and t — J models widely believed to model high- Tc superconductors. One of the universal features of high- Tc 
superconductors is the close proximity between the antiferromagnetic (AF) and the d— wave superconducting (dSC) 
phases [^-^. While an AF insulator appears to be diagonally opposite to a superconductor, their close proximity 
lead one of us (SCZ) to conjecture that they are in fact intimately related by an 5*0(5) symmetry group which unifies 
them 1^ . In this theory, the AF and dSC order parameters are grouped into a single five component vector Ua called 
superspin. The transition from AF to dSC is viewed as a superspin flop transition as a function of the chemical 
potential or doping, where the direction of the superspin changes abruptly. 

While the 5*0(5) symmetry was originally proposed in the context of an effective field-theory description of the high- 
Tc superconductors, its predictions can actually be tested within microscopic models. First numerical evidence for 
the approximate 50(5) symmetry of the Hubbard model came recently from exact diagonalizations of small-sized (10 
sites) clusters ||] studying dynamic correlations functions involving the AF/dSC rotation (tt) operators (l^,^. In this 
work we shall use the t — J model because of its more limited Hilbert space. Since the t — J model explicitly projects out 
the states in the upper Hubbard band, some of the questions ]TT[-p^ raised recently about the compatibility between 
the Mott Hubbard gap and 50(5) symmetry can also be answered explicitly. We use a general and direct recipe for 
checking microscopic Hamiltonians for 50(5) symmetry, i.e. the concept of "superspin multiplets". In particular, if 
there is an approximate 50(5) symmetry of the microscopic model, the low-energy states of this model should fall into 
irreps of 50(5). In a given quantum mechanical system, the direction of the 50(5) superspin vector is quantized in a 
way similar to an ordinary 50(3) spin, and the classically intuitive picture of the precession of the 50(5) superspin 
vector under the influence of the chemical potential [ppH can be identifled with the equal level-spacing between the 
members of 50(5) multiplets carrying different charge. Therefore, numerically identifying the low-lying states of the 
microscopic model with the 50(5) irreps can lead to detailed understanding of the one-to-one correspondence and 
the level crossing between the excited states of the AF and the dSC states, and thereby lead us to the microscopic 
mechanism by which the AF changes into the dSC state. While flnite size-calculations can not generally be used to 
prove the existence of long range order in infinite systems, the spectroscopic information about the 50(5) symmetry 
can be used as input for the effective field theory |p|, p5yi^ ] which captures the low energy and long distance physics 
of the problem. 

SO(5) Superspin Multipets, a Pyramid of Diamonds: The 50(5) Lie algebra is generated by ten operators 
Lab with a,b — 1, .., 5 and a < b. They obey the following commutation relation: 

[Lab, Led] — i{SacLbd + SbdLac — ^adLbc " SbcLad) (1) 

50(5) is a rank two algebra, we can therefore choose total the charge Q = Li^ ~ ^i-^e — M) and z component of the 
total spin Sz = —L23 to be the members of the Cartan subalgebra of mutually commuting generators. Here N^, the 
number of electrons and M, the number of lattice sites are both taken to be even. Moreover, the Casimir operator 
C = X)a<fc ^ab commutes with all generators and can be used to label the level of a representation. The tt operators 
are defined as follows 

TTa = Xl5(R)Cp-+Q(CTaCrj,)c_p-, (2) 

p 

where annihilates an electron with momentum p (we are suppressing the spin index) aa is the vector of Pauli 
matrices, and Q — (7r,7r) is the antiferromagnetic wave vector. If one takes g{p} — sgn[cospx — cospy), the 50(5) 
algebra closes exactly [T^ . However, for cluster calculations it is often more convenient to take g{p) — cospx — cospy, 
the numerical difference between these two choices is small. Together with the total spin raising and lowering operators 
5*, TTa and tt^ form the root generators of 50(5) and rotate different members of a multiplet into each other. 

In this paper, we are concerned with the tensorial representations of 50(5). Tensors with given symmetry types 
under permutation of their indices are classified by the Young tableaux p8[ . For the 50(5) group, tensors which 
have more than two antisymmetric indices can always be mapped to tensors with less or equal to two antisymmetric 
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indices by the invariant tensor e"''"^'^^, the fuhy antisymmetric index in 5 dimensions. Therefore, all tensorial irreps of 
SO{5) are characterized by two integers {u, v'), corresponding to the length of two rows in the Young tableaux [ p^ . 
The general {v, v') series can only be constructed from two different SO{h) vectors. However, as we shall see later, 
the low lying states of the t — J model can all be classified according to the restricted irreps [v^ 0) generated by the 
superspin vector alone. Therefore, we shall restrict ourselves to the fully symmetric tensors -F'ai,a2,..,a„ series (i^, 0) 
generated by the products of the 50(5) vector ria, satisfying 

[Lab, ric] = -iSbcUa + iSacm- 

Here Ua is the five-dimensional vector (A^ + A, S^, —i{A'' — A)), where A — (i/2) X^p'CpCr^C-p- denotes the dj;2_y2 
superconducting order parameter and S^^ — '^jjcL^J^Cp denotes the antiferromagnetic Neel vector . However, these 
representations are in general not irreducible. Since SO{N) transformations preserve the norm of a vector, the 
pairwise trace components of F should be projected out to obtain an irreducible tensor, i.e. i^ai,ai,..,a„ = 0. Since F 
is symmetric, the vanishing of the first pairwise trace ensures the vanishing for all other pairwise traces. Therefore, a 
pairwise traceless symmetric tensor has f ^ 1 \ f N + (i^ 2) 1\ ^^Qj^p^j^^j^^g^ which gives the dimension 
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of a level (i/, 0) (or simply v) irreps. For 5*0(5) we obtain = \{v + + 2){2v + 3), while for 50(3) 
reduces to the familiar degeneracy 2v + 1. The Casimir operator takes the value v{i' + 3) for a level v irreps. 

The linear combinations of tIq, A^ — rii ± in^, 5^ = n2 ± in^ and 5^, = are eigenvectors of Q and 5^, and 
can be used to label each component of a given irreps in the two dimensional coordinate space of Q and Sz- The 
diagrams of the multiplets take the form of a diamond as plotted in Fig. 1. Generally, a rank v irreps contains many 
spin multiplets, with total spin S ^ v multiplet being the largest member. At the top of the diamond, Q = v is, a 
spin singlet, at the next sublevel, Q — 1 is a spin triplet, the Q = v — 2 sublevel contains both a spin singlet and 
a quintet. Generally, the Q — v — p sublevel contains total spin 5 = p,p — 2,p — 4, .. multiplets. These different spin 
multiplets takes the form of nested diamonds in a multiplet. The different diamonds at level v are stacked together 
to form a pyramid, with v = Q singlet at the apex, and the v = M/2 diamond at the base of the pyramid. Each 
member of a given irreps is a box containing many microscopic states with the identical transformation properties 
under 50(5). 

AF and dSC states in the SO(5) supermultiplet: If we were dealing with a microscopic model with exact 
50(5) symmetry all the states at a given level v are degenerate with each other at any finite systems. Degeneracy 
between different v multiplets can only occur in a infinite system, signaling spontaneous symmetry breaking (SSB). 
On a finite system, the ground state is an 50(5) singlet, lying in the v = Q box at the apex of the 50(5) pyramid. The 
tendency towards SSB in the large system limit can be recognized from the scaling of the energies of excited states 
with higher irreps. For example, in the infinite size limit, a AF state with Neel vector in the xy plane is constructed 
from the linear superposition of states in the center column of the 50(5)-pyramid, while a dSC state is constructed 
from a linear superposition of states on the two (dSC) ridges of the 50(5)-pyramid. Because of the degeneracy within 
all multiplets, the AF and the dSC states constructed in the large system limit would have the same ground state 
energy. If one applies a chemical potential term i?^ , the members of a given multiplet with different charge quantum 
numbers will be linearly shifted by the ^2fiQ term, leading to an equal level difference within a given multiplet. 

In the microscopic Hubbard or the i — J model used to describe the high-Tc superconductors, the 50(5) symmetry 
is not exact, and there are different types of symmetry breaking terms. The symmetry breaking terms can also be 
classified according to irreducible tensors of the 50(5) Lie algebra. The chemical potential term mentioned above 
belongs to the 10-dimensional adjoint representation. The next simplest type of symmetry breaking term preserving 
spin rotation and charge conservation would be the Q = Sz = member of a 14 dimensional traceless and symmetric 
(2, 0) tensor, Hg, transforming like A+A^ — 2/3S^ ■ S^. This type of symmetry breaking has two important effects. 
First it can lead to mixing of states with ly quantum numbers differing by 2. The second more important effect is 
the removal of the degeneracy between the members of a supermultiplet carrying different charge quantum numbers. 
However, unlike the chemical potential term, it preserves the symmetry between the charge states with the same 
magnitude |Q|. This type of symmetry breaking term can remove the degeneracy between the AF and dSC states 
when fj, = 0, leading to a charge gap while keeping the spin excitations at low energy. However, with an applied 
chemical potential, the effects of these two types of symmetry breaking term can compensate each other for one type 
of charge states, say hole states with Q < 0, and there is a critical chemical potential /ic at which the multiplets with 
different charges Q < can recover their near degeneracy. As we shall see later, our overall numerical results can be 
consistently interpreted by these two types of explicit symmetry breaking terms. The competition between these two 
types of symmetry breaking terms is analogous to the competition between the spin anistropy and an applied uniform 
magnetic field in an antiferromagnet, as illustrated in reference H. 
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It is important to point out that, although Hg and can nearly compensate each other on one side of the charge 
states, the full 50(5) symmetry between low energy states of different signs of charge Q can not be recovered. In 
the t — J model for example, all states with Q > are projected out of the Hilbert space, and the 5*0(5) symmetry 
can only be approximately realized between the members of the supermultiplets on the Q < half of the 50(5) 
pyramid. However, Q < Q states are the relevant low-energy degrees of freedom in question, and the approximate 
50(5) symmetry between these states is sufficient to understand the full effect of doping. In this formalism we see 
the fundamental importance of the Mott-Hubbard gap (projecting out the Q > states) on the interplay between AF 
and dSC. Assuming that all the superspin multiplets are degenerate for the Q <Q states at a given /Xc, one can either 
form a AF ordered state by the linear superposition of the Q = members of the different v multiplets, or use the 
same coefficients to form a pure dSC ordered state by the linear superposition of the Q — —v members of the different 
V multiplets. These two states and the intermediate states between them Q are degenerate at [ic- However, since a 
macroscopic number of Q < states are used to form a phase-coherent pure dSC state, it has a finite hole density. 
Therefore, Mott insulating behavior at half-filling is compatible with the 50(5) symmetry: AF and dSC states are 
nearly degenerate at /ic, but they have different density. 

Results of exact diagonalization of the t-J model: Numerical demonstration of SSB requires careful study 
of the level spacing as a function of system sizes, and true SSB is generally hard to establish. However, it is relatively 
easy to recognize a nearly degenerate multiplet structure on a finite system. For example, if we have a weakly 
anisotropic Heisenberg model on a lattice, the energy splitting within a multiplet would be small compared to the 
splitting between multiplets. Although both level spacings may scale as 1/Af in the limit of large system size, the ratio 
of their difference can be independent of the system size and its smallness thus can be recognized even with limited 
finite-size data. In the following, we study the t — J model, the simplest model Hamiltonian which incorporates the 
key features of the strong-correlation limit: 

where denotes a summation over all nearest neighbor pairs on a 2D square lattice, V projects onto the subspace 
of states with no doubly occupied sites. The latter constraint reflects the strong correlations in the t//t — >• oo limit 
of the Hubbard model. The parameters —t and J are the nearest neighbor hopping and exchange integral. We have 
numerically diagonalized the t — J model on finite lattices of 16 and 18 sites (see Ref. Q for pictorial representations of 
these standard systems) and studied its low lying eigenstates up to total spin 3 and 6 holes. In addition to their spin 
and charge quantum numbers, these states are also labeled by their total momentum and the point group symmetry. 
In figure |l| we show how some of the low lying states of the 18-site cluster t—J model with J/t — 0.5 fit into the irreps 
of 50(5), up to the v = 3 supermultiplets. The ground states within the respective hole-number sector are labeled by 
an asterisk (we note that up to now this assignment of multiplets is only a conjecture; below we will present numerical 
evidence that these groups of states indeed are 50(5) multiplets). We see that all the different quantum numbers of 
the states are naturally accounted for by the quantum numbers of the superspin, and the levels with different charge 
Q are approximately equally spaced. More precisely, the mean level spacing within each multiplet up to Q = —2 is 
—2.9886 with a standard deviation of 0.0769. Therefore, at a chemical potential comparable to the mean level spacing, 
the superspin multiplets are nearly degenerate. There are four states inside the "nested diamonds" which cannot yet 
be fully identified with our current methods, they are marked by the () symbol. 

Now we wish to demonstrate that the different states inside a given multiplet can indeed be rotated into each other 
by the 50(5) root generators. In particular, we would like to show explicitly how higher spin AF states are rotated 
into the dSC states. Let us first briefly discuss the selection rules within flgure |l| The tt and tt^ operators act as 
raising and lowering operators and obey the selection rule Ai' — 0. In other words, they give transitions within the 
'diamonds' in figure |l|. In the presence of an Hg-type of perturbation Av — 2 transitions are also possible, but are 
expected to have smaller amplitude. Next, as mentioned above, the 5 operators A, and together form an 
50(5) vector. Consequently they play a role analogous to the dipole operator in 50(3). They obey the selection rule 
Ai/ = ±1, i.e. they allow transitions between the different diamonds. Thereby A^ is a spin singlet, transfers zero 
momentum and has Bi symmetry, whereas 5^ is a spin triplet, transfers momentum (tt, tt) and has Ai symmetry. 

is the operator relevant for neutron scattering, this experiment thus probes Aiy ~ ±1 transitions. 

As a 'diagnostic tool' to judge if a transition from a given state |^) by the operator O (which can be tt, Tr^, A or 5^) 
is possible, we compute the spectral function 

A^iu;) = 3i(*|6t I d\^) (3) 

TT U) — [H — i^ref) ^ *0+ 
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where 3 denotes the imaginary part and Eref is a suitably chosen reference energy. For finite systems spectral 
functions of the type (||) can be calculated exactly by means of the Lanczos algorithm An intense and isolated 
low energy peak in (say) the 7r-spectrum (^) then indicates that there is a state with Q + 2 into which is being 
transformed by tt. This should hold at least as long as the doped state is sufficiently low in energy to be still 
within the range of validity of the approximate SO(b) symmetry. Our strategy for the following, therefore, is to 
work ourselves through Figure |l| and to check '5'0(5)-allowed' and 'S'0(5)-forbidden' transitions by computing the 
respective spectral functions. We can then also investigate the influence of perturbations which could possibly break 
the 5*0(5) symmetry by studying their influence on the spectra, to see e.g. how 'S'0(5)-forbiddcn' transitions are 
enhanced by the perturbation. 

We now begin to discuss the results of our 'computer spectroscopy'. Figure |^ compares the spin correlation function 
at half-filling to the spectra of the Trt-operator for the lowest 2-hole states with d-symmetry and total spin and 2. 
The reference energy Eref has been taken as the energy of the half-filled ground state. The spin correlation function 
has a single dominant low-energy peak at an excitation energy of ^ 0.5 J (marked by an arrow) which clearly should 
be associated with a magnon state (all cluster states are exact eigenstates of so there is no SSB in the small cluster 

- which explains the finite excitation energy). The spectrum of the 7r''^-operator for the 2 hole, ^Bi state also shows a 
single, high intensity peak, which coincides with that of the spin correlation function, i.e. the excitation energies agree 
within computer accuracy (10~^^). Obviously the final states are identical, which shows that the Trt-operator indeed 
produces the spin resonance. Next, the spectrum of the tt '''-operator for the 2 hole, ^Bi state has a strong peak at high 
energy, plus a low-energy peak with significantly lower intensity, which again coincides with the spin resonance. Here 
it should be noted that the transition from the ^Bi state with momentum (0,0) to the ^Ai state with momentum 
(tt, tt) has Av — 2, it is therefore an 5*0(5) forbidden transition, but it is allowed by spin, momentum or point group 
symmetry. Inspection shows, that the intense high energy peak in the ^Bi spectrum is nothing but the lowest Ai 
state - this Ai/ = transition obviously is allowed by the S*0(5) selection rule, see Figure 0. The 5*0(5) selection 
rule thus is obeyed approximately, with the ratio of the two peaks in the ^Bi spectrum being a rough measure for 
the degree of symmetry breaking. The pattern of the explicit symmetry breaking is consistent with that of a second 
rank 5*0(5) tensor Hg. The intensity of the peaks in the various V-spectra' decreases rapidly with decreasing J/t 

- this indicates that corrections to the 7r-operator become more important at smaller J/t. On the other hand the 
additional peaks at higher energy in the 7r-spectra decrease rapidly as well and always stay well-separated in energy - 
restricting the Hilbert space to states below a cut-off frequency « 2 J would therefore give a very good eigenoperator 
of the Hamiltonian. 

We proceed to the 2-hole subspace,i.e. Q = —I. Figure ^ shows the spin correlation function for two holes, as well 
as various spectra of the tt and tt^ operator, Eref is the energy of the two-hole ground state. To begin with, the spin 
correlation function again has a dominant low energy peak, whose excitation energy scales approximately with J. 
The final state responsible for this peak is the lowest ^Bi state with momentum (tt, tt). Then, the spectrum of the tt 
operator calculated for the undoped ^Ai, k = (0, 0) state, and the spectrum of the operator for the ^Aik = (0, 0) 
state with 4 holes also do have intense low energy peaks. These peaks are well separated from some incoherent 
high-energy continua, which start above a lower bound of ~ 2 J, and again coincide to computer accuracy with the 
'^Bi state observed in the spin correlation function. This again confirms the interpretation of the spin resonance 
as a 'TT-excitation'. Looking at Figure |l| it becomes obvious that these two transitions have Ai^ = 0, i.e. they are 
'5*0(5)-allowed'. On the other hand, the 7r-spectrum for the undoped ground state, ^Ai, has a weaker peak at the 
position of the spin resonance. This Av = 2 transition is forbidden by the ideal 5*0(5) symmetry (see Figure]^), 
indicating again a weak breaking of the 5*0(5) symmetry. The decrease of the 'ir-peaks' with decreasing J/t is quite 
analogous as in the case of half-filled final states. The only exception are the peaks in the tt''' spectra (i.e. with initial 
states in the 4-hole subspace), which have practically zero weight for smaller J/t = 0.25. 

This also becomes clear if we study spectra with final states in the Q — —2 subspace. Figure ^ shows the spin 
correlation function at 4 holes, together with spectra of the tt operator for the ^Bi and ^Bi states of two holes with 
momentum (0, 0). ^From Figure |^, we see that the transition from the singlet state is forbidden, that from the quintet 
is allowed. Then, looking at Figure |l| it is apparent that there occurs a drastic change for J/t smaller than a cluster 
dependent value. For J/t > 0.5 in the 18-site cluster {J/t > 1 in the 16-site cluster) we have the 'standard situation': 
the dominant low energy '7r-peak' for the ^Bi initial state is more intense than that for the ^Bi state, indicating again 
a weak breaking of the 50(5) symmetry. Both '7r-peaks' coincide with the dominant low energy peak in the spin 
correlation function, which in turn stems from the lowest ^Ai state at (tTjTt), which confirms that the interpretation 
of this peak as a '7r-resonance' is valid throughout the low doping regime. On the other hand, for smaller J/t the 
correspondence between the spin correlation function and the 7r-spectra is essentially lost, in the case of the 16-site 
cluster the spin correlation function does not have a distinguishable low energy peak at all. Quite obviously, we have 
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reached the hmit of apphcabihty of the SO{5) symmetry, which seems to occur at a doping level of « 0.25%, with 
some dependence on the ratio J/t. This is roughly the same parameter range where dSC correlations vanish on the 
finite size cluster. 

Summarizing the study of the spin correlation function, we may say that the data are in overall agreement with an 
approximate 5*0(5) symmetry, in that 'S'0(5)-allowed transitions' usually have a larger intensity than the 'forbidden' 
ones. The data also show that the dominant low energy spin excitation at (tt, tt) always can be generated by adding 
or removing two electrons from the system by means of the 7r-operator, which obviously supports the conjecture of 
Demler and Zhang ]lO[ | that this low energy resonance in the dynamical spin correlation function is the hallmark of 
the approximate 5*0(5) symmetry. The agreement with the S0{5) symmetry deteriorates for higher doping levels 
and/or smaller J/t. 

We now proceed to map some additional transitions within the 5*0(5) multiplets. Figure ^ shows the spectrum of the 
TT-operator for the lowest triplet state with 2 holes (this state is the one which gives rise to the prominent peak in the 
spin correlation function in Figure ||). The initial state thus is '^Bi{n, tt) with 2 holes, and on the basis of Figure we 
expect a strong transition to the ""^^1(0, 0) state with 4 holes (i.e. the 4 hole ground state). As the reference energy we 
choose the energy of the 4 hole ground state, and Figure H then clearly shows a pronounced peak in the 7r-spectrum 
with zero excitation energy, precisely as expected on the basis of the 5*0(5) symmetry. Figures ^-|^ thus demonstrate 
that the 4-hole ground state can be obtained by 2-fold '7r-rotation' from the lowest ^Ai state with momentum (0, 0) 
at half- filling. Similarly, the 2-hole ground state can be obtained by tt- rotation from the lowest half- filled ^Ai state 
with momentum (7r,7r). 

To summarize this section, we have shown that the transitions induced by the tt and S{Q) operators can be well 
understood in the framework of a weakly broken 5*0(5) symmetry. We have explicitly identified the multiplets with 
v — 0,1,2,3 and shown that the 7r-operator gives transitions between the members of a multiplet with different 
Q. Moreover, we have verified explicitly that the spin correlation function 'operates in the same subspace' as the tt- 
operator, and that the prominent low energy peaks in the dynamical spin correlation function corresponds to members 
of the 50(5) multiplets for all dopings < 25%. 

Conclusions: We have numerically diagonalized the low lying states of the t — J model near half-filling and 
found that they fit into irreps of 5*0(5) symmetry group. At a critical value of the chemical potential, the superspin 
multiplets are nearly degenerate, and therefore higher spin AF states at half-filling can be freely rotated into dSC 
states away from half-filling. There are clearly visible effects of 50(5) symmetry breaking, which to the lowest 
order can be identified with the type of a symmetric traceless rank two tensor. Our overall result suggest that the 
low-energy dynamics of the t — J model can be described by a quantum 50(5) nonlinear a model with anisotropic 
couplings, and the transition from AF to dSC phase can be identified with that of a superspin flop transition It 
is truely remarkable that while the physical properties of AF and dSC states are intrinsically different and they are 
characterized by very different forms of order, there exists nevertheless a fundamental 5*0(5) symmetry that unifies 
them. The dichotomy between their apparent difference and fundamental unity is in our view a key which can unlock 
the mystery of the high Tc superconductivity. 
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Appendix: Perturbations to the t~J Model and Their Influence on the Approximate 50(5) Symmetry 
Recently, Baskaran and Anderson ||l3|] raised some questions concerning the effect of the diagonal hopping and nearest 
neighbor Coulomb interaction on the approximate 5*0(5) symmetry. It is then of important to check whether these 
perturbations are essentially irrelevant or they could lead to a breakdown of the (approximate) 5*0(5) symmetry. 
Again, we resort to exact diagonalization calculations to address this question. To begin with, we consider the effect of 
an extra 'Coulomb repulsion' between holes on nearest neighbors. More precisely, we add the term Hy = ^ j) ''^i^j 
to the Hamiltonian, where denotes the electron density on site i. Figure |^ then compares some spectra of the tt-, 
TT^-operator and the zero momentum pair-operator A for different values of V. The left panel shows spectra with final 
states in the two-hole sector, the reference energy is that of the half-filled ground state. With this choice the excitation 
energies of the dominant low-energy peaks in both spectra increase with V. This is natural because the repulsion 
V is not operative at half-filling, but will tend to increase the energies of hole-doped states. The increase, however, 
is significantly less than expected, being only approximately 0.5t for V—2t. This can hardly come as a surprise, 
because we have dEo/dV = A{nini+±)^ i.e. the nearest neighbor density correlation function of holes. The latter 
quantity is quite small in the physical range of parameters, so V does not have a great impact. More importantly, 
the difference of the excitation energies of the A and tt operators is practically independent of V . This difference of 
excitation energies would give the energy required to remove a fc = (0, 0), d-wave singlet pair from the system and 
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reinsert a fc = (tt, tt) d-wave triplet pair. This is exactly what is done in a neutron-scattering experiment - a Cooper 
pair from the condensate is turned into a 7r-pair while scattering an incoming neutron. The energy difference of the 
peak energies in Figure ^ thus should correspond to the energy of the peak in the inelastic neutron scattering cross 
section, and Figure |^ clearly shows that even rather strong repulsion between the holes leaves this energy unchanged. 
Moreover, we note that the weight of the peaks decreases only slightly with V - the decrease is also very similar for 
the A and tt operator; this would suggest that as long as superconductivity 'survives' the influence of V, so does 
the TT-resonance. It should also be noted that the figure actually compares an 'iS'0(5) forbidden' transition (from the 
half-filled ^Ai{0,0) state to the ^i3i(7r, tt) with two holes), and an 'SO{5) allowed' transition (from the ^Ai(0, 0) 4 
hole ground state to the ^i3i(7r, tt) state). The ratio of intensities for both transitions is !=s 1/3 and it stays so more 
or less independently of V. This indicates that the degree of symmetry breaking is not affected significantly by V. 
Next, the right-hand panel in Figure ^ compares the tt and A spectra calculated for the ground state with two holes 
- it shows similar features, in particular the difference of excitation energies is independent of V, and the weights of 
the peaks decrease in a very similar fashion with V. 

We now consider the influence of a next-nearest neighbor hopping integral t'. We choose a t' between (l,l)-like 
neighbors with opposite sign as t; for noninteracting electrons this would produce a Fermi surface similar to the LDA 
predictions. Figure again compares the spectra of the 7r-operator and the A-operator with different t' and different 
doping levels. The reference energy again is the ground state energy at half-filling in the left panel, which shows 
spectra with final states in the Q — —1 subspace; in the right panel, which shows final states in the Q — —2 subspace 
Eref is the ground state energy of two holes. The overall picture is comparable to that seen in Figure ^, i.e. the 
difference in excitation energies is nearly independent of t' /t, and in fact even decreases with increasing t' /t. In the 
spectra with 4-hole final states this obviously leads even to a kind of level-crossing, in that the lowest '7r-peak' comes 
down below the lowest A-peak for large t' . The intensity of both low-energy peaks decreases in an essentially similar 
fashion with increasing \t' /t\. One can, however, realize a kind of 'crossover' between \t' /t\ = 0.1 and \t' /t\ = 0.2, 
where the spectral weight of the 7r^-spectra drops sharply. The ultimate reason is a level crossing in the 4 hole sector 
from the 'S'0(5)-compatible' ^Ai ground state to a ^Bi ground state, which occurs in between these two values of 
\t' /t\. We note that 4 holes in 18 sites correspond to a hole density of 22%, which is nominally far overdoped - the 
drop in the Tr^-spectra thus is not really a reason for concern. Moreover, the peak in the spectrum of the 7r-operator 
(which is not affected by the level crossing) stays well defined, its intensity decreasing slightly and in proportion to 
that of the 'A-peak'. All in all it is obvious that larger values of t' degrade the 5*0(5) symmetry. On the other hand, 
practically all of our data show an intimate relationship between the d-wave pairing amplitude and the 7r-amplitude. 
If the TT resonance is suppressed, be it due to high doping, large V or large t\ so is usually the d-wave pairing. It is 
then only natural to conclude that as long as the d-wave pairing 'survives' the influence of perturbations, so does the 
TT-resonance. 
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FIG. 1. The upper diagram illustrates a general level v irreps of 5*0(5). Every state can be labeled by Q and S^- The 
maximal charge is Q = ±!/. The states labeled by a x form the shape of a diamond, while states inside the nested diamonds 
are labeled by o and A. Overlapping states with same Q and Sz are distinguished by their S quantum numbers. The lower 
diagrams are for v = 1,2,3 irreps of SO{5). The figure shows the energies of some low energy states for the 18-site cluster 
with J/t — 0.5. The states are grouped into different multiplets and are labeled by the spin, point group symmetry, and total 
momentum. Ai denotes the totally symmetric, Bi the ci^2_y2-like representation of the C4„ symmetry group. The () symbol 
denotes as yet unidentified members of the respective multiplet. 

FIG. 2. Comparison of spectral functions with undoped final states: dynamical spin correlation function for momentum 
transfer Q, calculated for the half-filled ^j4i(0, 0) ground state; spectrum of the Tr^-operator, calculated for the ^Bi{0, 0) ground 
state in the Q — —1 sector; spectrum of the vr^-operator, calculated for the lowest ^Bi{0, 0) state with Q = —1. Data are shown 
for different cluster sizes and values of the ratio J/t. 

FIG. 3. Spectral functions with final states in the Q = — 1 subspace: dynamical spin correlation function for momentum 
transfer Q, calculated for the ^i3i(0,0) ground state; spectrum of the 7r^-operator, calculated for the ^Ai(0,0) ground state 
in the Q = —2 sector; spectra of the vr-operator, calculated for the half-filled ^Ai{0,0) ground state and the lowest half-filled 
^Ai{0,0) state. 

FIG. 4. Spectral functions with final states in the Q = —2 subspace: dynamical spin correlation function for momentum 
transfer Q, calculated for the ^^i(0, 0) ground state; spectra of the vr-operator, calculated for the '^i3i(0, 0) ground state and 
the lowest ^Bi{0, 0) state in the Q — ~l subspace. 

FIG. 5. Spectrum of the vr-operator in the 'spin resonance state' ^i3i(7r,7r) at Q = —1; this state is the final state corre- 
sponding to the dominant peak in Figure 0. 



FIG. 6. Spectra of the n and A-operator for different V. 



FIG. 7. Spectra of the tt and A-operator for different \t' /t\. t' has opposite sign as t. 
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